Non-extensive statistics and the stellar polytrope index 
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Polytropic distribution functions provide the simplest 
models for self-gravitating stellar systems These 
models are characterized by an equation of the state of 
the form p = Kp'^ , where K and 7 are constants. The 
polytropic index n is defined by the following expression: 
7 = 1 -|- i which in the limit n — > 00 or, equivalently, 
7 = 1 results in an equation of state identical to the one 
of an isothermal body of gases, i.e., p — Kp. The distri- 
bution function associated to these systems is given by 
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where pi and a — {ksT /mY^'^ are, respectively, the den- 
sity and the velocity dispersion. The quantity e is equal 
to a constant minus the total energy per unit mass of 
an individual star moving in the global galactic gravita- 
tional potential. In fact, e is the relative energy of a star 
defined by e = — /2, where ^ is the relative poten- 
tial, given by \1/ = 00 — <^ is the star potential energy 
and is a arbitrary constant. In particular, the inte- 
grability condition J /d^v < 00 results in the constraint 
n > 1/2 (for a complete study on the gaseous polytropes 
and isothermal gas spheres, see |3|)- 

On the other hand, a non-extensive statistical formal- 
ism has been proposed as a possible extension of the clas- 
sical one, which intends to study more properly systems 
that possess long-range interactions. Such a formalism is 
based on the entropy formula 
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where pi is the probability of the ith microstate and q 
is the non-extensive parameter. The Botzmann-Gibbs 
extensive formula. Si = — fc^ J2iPi l^^Pii is readily recov- 
ered in the extensive limit q = I and the g-logarithmic 
function above is defined as 4] 

ln,/ = ^^ , />0. (3) 
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Some attempts to build a kinetic counterpart for this 
non-extensive statistical formalism has been recently pro- 
posed in Refs. 1^. The main result of these works 
consists in showing that this extended formulation leads 
to a new velocity distribution for free particles given by 
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where the quantity Bq is a ^-dependent normalization 
constant whose expressions for q < 1 and q > I are re- 
spectively given by 
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where the coefficients G and Ti. are given hy Q = np{l — 
g)i/2(^)(^) and H = np{q - 1)^/2, T is the tem- 
perature and Up is the particle number density. This 
q-distribution can be derived at least from two differ- 
ent methods, namely, a generalization of the Maxwell 
ansatz, f{v) ^ f{vx)f{vy)f{vz), which follows from the 
introduction of statistical correlations between the com- 
ponents of the velocities ^ and a new formulation for 
the Boltzmann i?-theorem, which requires g > 0]. 

In the astrophysical domain, both the non-extensive 
statistical formalism and its kinetic counterpart have 
been applied in a wide range of problems. The very first 
application of this non-extensive statistics was done in 
connection with stellar polytropes j^, with several au- 
thors suggesting different expressions for the polytropic 
index n as a function of the non-extensive parameter q 
IE M> ^1 ^1 ■ Recently, a numerical study of stellar dy- 
namical evolution for self-gravitating systems was per- 
formed by Taruya and Sakagami 13] while Chavanis 
considered the Tsallis' entropy as a particular iJ-function 
corresponding to isothermal stellar systems and stellar 
polytropes. In this latter approach, the maximization 
of the _ff-function at fixed mass and energy reveals a 
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thermodynamical analogy with the study of the dynam- 
ical stability in coUisionless stellar systems. The non- 
extensive kinetic approach has been used to study the 
Jeans' gravitational instability [T3.ll5l |. 

Another recent application of this non-extensive ki- 
netic theory was to investigate the negative nature of 
the heat capacity for isolated self-gravitating system . 
This study was performed by considering the follow- 
ing steps: first, consider a simple analogy between self- 
gravitating system and ideal gas (IG) or, in other words, 
the equivalence between IG kinetical and internal ener- 
gies, given by < >— ^ksT, where m is the mass 
of a particle (e.g., a star). If such a system is composed 
by N particles, its total kinetic energy is K = ^NksT 
which, according to the viral theorem, is equal to minus 
total energy, i.e., E = —K. Therefore, the heat capacity 
of the system is given by (see for details) 
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Second, in order to investigate the g-dependence of the 
heat capacity, we consider a cloud of ideal gas within the 
non-relativistic gravitational context, which is analog to 
a self-gravitating coUisionless gas. The kinetic energy of 
this system is simply given hy K — < > , with the 
statistical content included in the average square velocity 
of the particles < >. Indeed, the non-cxtcnsivity may 
be introduced through a new derivation of expectation 
value j21|| 
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FIG. 1: The quantity Cv/Nks for a self-gravitating collision- 
less gas obeying the non-extensive Tsallis' g-statistic is shown 
as a function of the polytropic index n. The polytropic index 
and non-extensive parameter q are related by Eq. 11211 . As in- 
dicated, the shadowed regions stand for the constraints from 
the second law of thermodynamics and the negative nature of 
Cv. 
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is a thermal cutoff on the 



maximum value allowed for the velocity of the particles 
(g < 1), whereas for the power law without cutoff {q > 1) 
Vra —>■ oo. This q-cxpectaiou value can be easily evalu- 
ated for q ^ I, resulting in 
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for q < 5/3. 
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Now, combining the g-expectation value for the square 
velocity of the particles with the definition of the heat 
capacity Cy one finds 
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which clearly places an upper limit to the non-extensive 
parameter, i.e., q < 5/3. Following this reasoning, it 
is therefore natural that our kinetic approach uses the 
free particle distribution, with the long range nature of 
gravity being introduced via the virial theorem (see 
for details). We note that, in the domain of ensemble 
theory, there are some analyses in the literature dealing 
with the classical IG within this nonextensive scenario 
[HEIIIIIS- In particular, Ref. [3 showed that the 



IG naturally exhibits negative heat capacity and that its 
pressure and particle number density obey a polytrope- 
type relation. Here, however, we assume that all IG prop- 
erties are accounted for by the g-Maxwellian (see 0, llg 
for details). 

In this Letter, by considering the different expressions 
relating the polytropic index n and the non-extensive pa- 
rameter q we discuss some possible physical constraints 
on n which arise from the second law of thermodynam- 
ics {q > 0) and the negativeness of the heat capacity 
{q < 5/3). We show that the existent expressions be- 
tween n and q (see H, |^ 0) can be naturally linked 
with the heat capacity imposing limits on the polytropic 
index. It is worth mentioning that the constraint q > 0, 
arising from the generalized 77-theorem, does not provide 
a strong case for limiting the polytropic index n. This 
corresponds to a weak condition of the generalized H- 
theorem since it is possible to show that there are two 
arrow of time for g > and g < 0. 

In order to investigate such constraints we first present 
the different expressions for the polytropic index n as a 
function of the non-extensive parameter g discussed in 
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C/NK =(9-6n)/4n 



TABLE I: Limits to n 



Thermodynamic constraints 



case (i) 



> 

< 5/3 



< n < i 



the polytropic index: 



and 



^ < n < I 
< n < I 



(ii) = ^^-^NkB. (13) 
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FIG. 2: The same as in Fig. 1 for the case (ii). Here the 
polytropic index and non-extensive parameter q are related 
byEq. CT . 



the literature, i.e., 



(i) ^=o + l [8]' (10) 

and 

(ii) "^2ir^ t^'^'^l' (^^^ 

The case (i) arises from an extension of Padmanab- 
han's classical analysis of the Antonov instability for the 
case of polytropic distribution function . By using ki- 
netic arguments of the Tsallis' non-extensive formalism, 
the case (ii) is obtained from an extension of the poly- 
tropic Lane-Emden spheres as well as from a study of 
maximum entropy solutions of the Vlasov-Poisson equa- 
tions describing self-gravitating systems. We emphasize 
that the polytropic family of solutions to the coUisionless 
Boltzmann equation (parameterized by the index n) is 
always the same, which means that expressions (jlOf) and 
constitute only different ways of parameterizing (in 
terms of non-extensive parameter q) this polytropic fam- 
ily. In both cases, Eqs. (|10|l and Hll|l . the same choice of 
the statistical average, the so-called g-expectation value, 
has been used. As expected, for all the above expres- 
sions the Maxwellian isothermal spheres are obtained in 
the limit q = 1 or, equivalently, n — oo. 

By combining Eqs. l(T(Hl - l(TT|l with Eq. Q, one obtains 
the following expressions relating the heat capacity and 



From the above expressions, it is straightforward to see 
that by imposing the constraints on the non-extensive pa- 
rameter q from the second law of thermodynamics {q > 0) 
and the negative nature of Cy (g < 5/3) it is possible to 
limit regions in the — n plane. To better visualize such 
constraints, in Figs. 1 and 2 we show, respectively for the 
cases (i) and (ii), the dimensionless quantity Cy/Nks as 
a function of the polytropic index n. 

In Figures 1 [case (i)] and 2 [case (ii)] we show that the 
constraint q > results, respectively, in the forbidden in- 
tervals for the polytropic index i < n < | and | < n < | 
while the limit resulting from Eq. jnj, i.e., <? < § implies 
< n < i and < n < |. In all cases the bound on q 
from the second law of thermodynamics (g > 0) restricts 
only negative regions in the plane Cy /Nkg — n while, as 
expected, the themodynamical limit C'v < forbids any 
possible positive region in the plots. In particular, we 
also note that the constraints on n from Eq. H12|l (case 
i) seem to be in agreement with the theorem proved by 
Antonov and Dremus Q who showed that the equilib- 
rium distribution is dynamically stable for values of the 
polytrope index n > 3/2 or, equivalently, g > (see Fig. 
1). The bounds on the polytopic index n obtained in this 
Letter are summarized in Table I. 

Finally, we emphasize that in a more quantitative case, 
i.e., in which a self-gravitating ideal gas is contained in 
a spherical container of radius R our constraints on the 
polytropic index n must be valid only in the negative re- 
gion of the heat capacity. As widely known, when the 
potential energy is taking into account or, equivalently, 
the virial theorem is modified to introduce a pressure 
term, 2K + W ^ p{R), the heat capacity can assume 
both positive and negative values. We also note that 
a possible connection with stability of spherical systems 
can be introduced through the Doremus-Feix-Baumann 
theorem which implies, for the extensive case, that poly- 
trojp es with ^ < n < | are unrealistic (see, however. 
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